Multivariate Wavelet Frames ^ 

M. Skopina 



Abstract 

We proved that for any matrix dilation and for any positive in- 
teger n, there exists a compactly supported tight wavelet frame with 
approximation order n. Explicit methods for construction of dual and 
tight wavelet frames with a given number of vanishing moments are 
suggested. 

1. Introduction 

Approximation order of a wavelet frame decomposition 

oo r 

E EE</.'^^>v$' (1) 

j=-oo fegZ'* 1^=1 

is closely related to the number of vanishing moments. It follows from Propo- 
sition 5.7 of [3] that a MRA-based compactly supported dual wavelet frame 
system with diagonal matrix dilation provides approximation order N when- 
ever the Fourier transform of all its generators has zero up to order N at 
the origin (in our terminology, this means that vanishing moment property 
of order — 1 holds for all the generators). We give a proof of another 
fact: for a wide class of dual wavelet frames ( not necessary MRA-based, 
not necessary with diagonal matrix dilation and not necessary compactly 
supported), vanishing moment property for ip^'^\ z/ = 1, . . . , r, is sufficient for 
the corresponding approximation order of (1). On the other hand, vanishing 
moment property of order for ip^'^\ v = 1, . . . ,r, is a necessary condition 
for the system {ipp}} to be a frame. This explains why vanishing moments 
property is useful for frames, and why study of this topic is important. Some 
additional arguments are presented in Section 2. 
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If a MRA-based wavelet system is a basis, then the number of its vanishing 
moments depends only on the dual generating refinable function. Situation is 
essentially different for frames. Two pairs of dual wavelet frames may be gen- 
erated by the same refinable functions and have different number of vanishing 
moments. The goal of this paper is to describe refinable functions generat- 
ing dual wavelet systems (potential frames) with vanishing moments and to 
present an explicit method for construction of compactly supported wavelet 
frames with arbitrary number of vanishing moments. A close problem was 
studied by Ming-Jun Lai and A. Petukhov[10] for univariate dual wavelet 
frames. Their technique is not appropriate for multi-dimensional investiga- 
tions because zero properties of multivariate masks can not be described by 
means of factorization in contrast to the one-dimensional case. 

Throughout the paper we will use the following notations. 

N is the set of positive integers, denotes the (i-dimensional Euclidean 
space, X = {xi, . . . , Xd), y = {yi, . . . ,yd) are its elements (vectors), {x, y) = 
Xiyi + - ■ ■+Xdyd, \x\ — Vi^c^, = (0, . . . , 1, . . . , 0) is the j-th unit vector in 
R*^, = (0, . . . , 0) e M'^, 1 = (1, . . . , 1) e K<^; Z*^ is the integer lattice in M'^. 
For x,y E W^, we write x > y iixj > yj, j — 1, . . . , d; Z'^ = {x e : x > 0}. 

Ua,f3 e Zl, a,b e R", we set a\ ^ fl a,!, h) = a' ^ fl a/^ 

j=i \ / j=i 

[a] — cvj, D^f — gaixi d^dx ' denotes Kronecker delta; T*^ is the unit 

d-dimensional torus; C is the set of complex numbers. 

Let M be a non-degenerate d x d integer matrix whose eigenvalues are 
bigger than 1 in module, M* is the conjugate matrix to M, denotes the unit 
dxd matrix. We say that numbers k,n & Z''' are congruent modulo M (write 
k = n (mod M)) ii k-n ^ Mi, £ e Z'^. The integer lattice Z'^ is spht into 
cosets with respect to the introduced relation of congruence. The number 
of cosets is equal to |detM| (see, e.g., [12, §2.2]). Let us take an arbitrary 
representative from each coset, call them digits and denote the set of digits by 
D[M) . Throughout the paper we consider that such a matrix M is fixed, m = 
I det M| , D{M) = {so, ■ ■ ■ , Sm-i}, so = 0, R{M) = {M'^o, . . . , M'^Sm-i}. 

Let ip^''\ip^''^ e L2(W^), u-l,...,r. Dual wavelet systems {V']?}, {V']?}, 
where := ''^i^^^^M^ ■ +k), ^^"^ := m^'/2^M(MJ- . +^)^ j^j. ^ X'^, 
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V — 1, . . . , r, form dual frames if for each / G L2 



jk 



(2) 



j=-oo 



and 



j = -0O A;eZ<* I/=l 

00 r 

^ii/ir< E EEK/'^]?)i'<^ii/f' ^'^>o- (4) 



Let a MRA in L2{W^) be generated by a scaling function (p which satisfies 
the refinement equation 

(p{x) = mo(M*-^x)lp{M*-^x), 

where mo G L2{T'^) is its mask (rcfinablc mask). For any rrii, G L2{T'^), 
there exists a unique set of functions /ij,fc G L2(T'^), k = 0, . . . ,r, (polyphase 
representatives of m,^) so that 

^ m— 1 

m,(x) = ^ E e'"'(^'='"V.ik(M*x). (5) 

fc=0 



The functions /Xj^fe can be expressed by 



It is clear from these formulas that a function rrii, is differentiable {n times) on 
R{M*) if and only if its polyphase representatives fi^k, /c = 0, . . . , m — 1, are 
differentiable (n times) at the origin and m,, is a trigonometric polynomial if 
and only if its polyphase representatives, are a trigonometric polynomials. 

Let now another MRA be generated by a scaling function Ip with a mask 
rho. According to Unitary Extension Principle [16], to construct dual wavelet 
frames one finds wavelet masks m,^, mj,, = 1, . . . , r, r > m — 1, so that the 
polyphase matrices 



M :-- 



A*0,m-1 
f^r,m—l J 



( 



M := 



A*oo 




satisfy 

M^M^Im, (6) 
and define wavelet functions by 

= m^{M*-^x)ip{M*-^x), 

The corresponding dual wavelet systems are {i^j'k}-> {'^jk} ^^^'^ 
generated by </?, (p (or MRA-based) . 

It is known that if = then {V'jfc^} is a tight frame in L2(M.'^). If 
A4,M are arbitrary matrixes satisfying (6), under some additional assump- 
tions on (p, ip,mv,mi, (see [6], [7], [12, §2.7]), we can state that {V'jfe^}, {V'jfe^} 
are dual frames in L2(]R'^). 

Definition 1 We say that a wavelet system {ipj'^'} has vanishing moments 

up to order a, a G Z'^, ( has VM^ property in the sequel), if D^ip^^^Q) = 0, 
u^l,...,r, for alip e 1.%, P<a. 

Usually it is more useful to control univariate order of vanishing moment 
property (for example, to apply Taylor formula). 

Definition 2 We say that a wavelet system {'4'jk} has vanishing moments 

up to order n, n e Z+, (has property in the sequel) if D^ip'^''\0) = 0, 

u^l,...,r, for alip e \j3] < n. 

2. Why vanishing moments are needed for frames? 

Theorem 3 Any wavelet frame {ij^jk} has VM'^ property provided all the 
functions jp^'^^ are bounded and continuous at the origin. 

Proof . Let {ipp}} be a wavelet frame in L2(M'^). Assume that ip^^^O) ^ 
for at least one i>. Take a function / G L2(M°') whose Fourier transform is 
compactly supported. Using standard technique, we have 

EK/'^S^)I' = / li^^'KM-^xl'lM'dx + Rj, 
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where 



Rj< f Yl + li^^^HM-^xl lip^^^M-^x + l)\dx. 



If j is large enough positive integer, then the supports of the functions / and 
/(• + MH) are disjoint for all leZ'^Jj^O. Hence 



This contradict to (3).0 

Let {i/jj'k}-, {'^jk} be dual frames. Decompositions (2) is said to have 
approximation order n if there exist C > and A > 1 such that for any 
function / in the Sobolev space W2 



< c 



WX' 



(7) 



Since all spectrum of the operator M"^ is located in the circle \x\ < r(M~^), 
where r{M~^) :— limj_^oo is the spectral radius of M~^, and there 

exists at least one point of the spectrum on the boundary of the circle, (7) 
follows from 



with A which is less than the minimal module of an eigenvalue of M. For a 
diagonal matrix M = cld, (7) holds with A = c. 

Theorem 4 Let {V'j-fe^}; {V'jfc^} dual wavelet frames, 

(1 + Fir 



for all v = 1, . . . ,r, and almost all x e M*^. // {ipj'^i)} has VM"' ^ property, 
then decompositions (2) have approximation order n. 

We need the following auxiliary statement for the proof. 
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Lemma 5 Let r] be a positive bounded function decreasing on [0, oo) so that 
rj{\x\) is summable on W^. Then there exists a constant K depending on rj 
and d such that 



^ r]{\x + k\)r]{\y + k\) < Kr] 



\x-y\ 



(8) 



8 

for all x,y e W'-. 

Proof . First of all note that there exist constants Ki, K2 such that 

v{\t\) < i^i, Yl ^(1^ + ^1) - ^2 Vt e M'^. 

Since both the left and the right hand sides of (8) are invariant with respect 
to the operation {x,y) — >• (x + l,y + I), I G Z'^, we can assume that < 
If \y\ < 2Vd, then 

5^ r]{\x + k\)r]{\y + k\) < K.K^ < - y\). 

Now let \y\ > 2\/d. Since jo; — |/| < 2\y\ and \y\ > 4:\x\, we have 

\x-y 



|fc|<M 



|fc|<M 



K2V 



J2 V{\x + k\)7]{\y + k\)< ^ ( in j ^(1^ + < i^^^r; I ^ 1 < 



|fc|>M 



|fc|>M 



\k\ 



K2V ( 1 -0 



Proof of Theorem 4. Let / G W^, j e Z+. It follows from (2) that 



2 »>j ^^=1 



(9) 
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Since VM'^ property is equivalent to 



using Taylor formula, we have 



f{y)4l\y)dy 



< 



1 

J ^^1^- ^11^1^(^)1 / \D''f{x + t{y-x))\dt. 



[a\ — n 



From this, due to Lemma 5 and Cauchy-Bunyakovskii inequality, we obtain 



< 



Cl / dx 



dy\x-y\-Y.\^'£^y)'l^^^^)\ I E + 



< 



C3 J dx J dy J dt 



{l + \M%x-y)\y 



(iw / dt 



Rd 



(1+ |M%|)' 



J2 \Dy{x + t(y-x))\'- 

\ [a]=n 

/ 

J2 \Dy{x + tu)\' 

\ [a]=n 



\ 



I I 



< 



/ dx 



du / dt 



1 

j du j dt— 



+ \M-'u\)^ 



mm 



(1 + \M'u\)^- 



J2 \D-f{x + tu)\' 



W2 





u 




(1 + 


|M*m|)^ 



< 



3 I \\J \\W^ 



\M- 





V 


n 


dv 


(1+ 


V 


)7 



< 



Cs I ||/||w2||M-"« 



1 =C4||/||^2||M-^ll2n 

(1 + I f 



It remains to combine this estimation with (9). <0 

The scheme of the proof of Theorem 4 can be apphed for some other 
approximation problems. In a similar way, it is possible to estimate coeffi- 
cients of decomposition (2), to find the order of approximation at a point for 
functions with some special local properties, e.g., for the class (introduced 
by Calderon and Zygmund [1]) of functions / e L(R'') such that. 



j \f{x)-P{x-xo)\dx^o{h''), 



\x—xo\<h 

where P is an algebraic polynomial. Detailed consideration of these problems 

is out of our today's interest. Our main goal is investigation of KM" property 
for frames and development of methods for their construction. 



3. Polyphase characterization of vanishing moments 
property 

Now we will consider only wavelet systems {tpjk}-) {i^jk} which are generated 
by scaling functions ip, (p whose masks mo, mo are continuous at the origin 
and mo(0) = mo(0) = 1. 

Assume that the functions ^, mi, . . . , m^ have derivatives up to order a 
at the origin. It easily follows from Leibniz formula that VMa property holds 
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if and only if 



D^{m^{M*-^x) 



x=0 



0,i/=l,...,r, V/?eZt/3< 



a. 



(10) 



In the case r — m — 1, there exist different criterions for vanishing moment. 
It is known [9] how to describe vanishing moment property in terms of hnear 

identities for Fourier coefficients of the dual refinable mask (so-called sum 
rule). Some other descriptions of masks providing VMa property are found 
in terms of zero-conditions [9] and in terms of containment in a quotient 
ideal [11]. The following polyphase criterion was given in [17]: VM^ property 



is valid for } if and only if there exist complex numbers A^, 7 e 
7 < a, such that Aq = 1, 



i^^/^ofe(0) = 




0<7</3 



-2TTiM-^Skf-"' V/? e Z^,/? < a, (11) 



for each k = 0,...,m — 1. The set of parameters A-y in (11) is unique, and 
A^ does not depend on a due to the following statement. 

Proposition 6 [17] If (11) is valid for the polyphase representatives of rfiQ, 
then 

(12) 



A^ = {mo{M*-^x)) 



x=0 



for all PelA, I3<a. 



So, in the case r 



m — 1, V Ma property for {'(/'j^.^} depends only on mo, 

i.e. only the first raw of the matrix M. is responsible for vanishing moments 
of wavelets generated by the matrix A^. In the case r > m—1, VM^ property 
for {V'jfe^} depends also on the way of construction of matrixes M, M. This 
may be illustrated by the following example. 



Let d^l, M 



m 



— 1 



A*oo — A*oi — A*oo — A*oi = 



1 

72 



1 

72 



M' 



'2 

V2 





73 / 



M' 



V 



1 

V2 
1 

72 
1 

72 



1 

1 

72 

72 / 



Either of pairs Ai.Ai and Ai', Ai' satisfies (6). The matrixes Ai,Ai generate 



wavelet masks mAx) — m2{x) — mi(x) — 7712(2;) 



1 

272 



272^ 



2mx 



It is clear 
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that mi(0) = m2(0) = mi(0) = m2(0) = 0, i.e. for the corresponding wavelet 

systems VMq property is vaUd. The matrixes Ai', A4' generate wavelet masks 
'm[{x) = |, m2{x) = ie^''*'', m'i(x) = | — ie^™,m2(x) = — | + |e^™, and 
we have m[{0) ^ 0,m'2(0) 7^ 0. 

Theorem 7 Le^ a G Z'^, r > m — 1, i/ie functions f^u,k,]iuk G -^2(T'^), 
z/, /c = 0, . . . ,r, have derivatives up to order a at the origin, the matrixes 
M ■■= {l^i^kK,k=o andN := {Jiyk}l,k=Q satisfy 

MW=Ir+i; (13) 

and let masks mQ,mi, . . . ,mm~i be defined by (5). Then condition (10) is 

valid if and only if 

(a) there exist G C, 7 G U^., 7 < a, such that Aq = 1 and (11) holds for 
k = 0, . . . ,m — 1; 

(b) D^Jlok{0) = 0, k = m, . . . ,r for all ^ e Z^, < a. 

Proof . Suppose that (10) is valid. We will prove (a) and (6) by induction 
on a. Check the initial step for a = 0. Let mi,{0) = 0, u = 1, . . . ,r. It follows 
form (5) that 

m— 1 

^/x,,(0)=0, u^l,...,r. (14) 

A:=0 

On the other hand, by (13), 

r 

X^/^ofc(0)/i^fc(0) = 0, 1/ = 1, . . . ,m - 1. 

fc=0 

Because of linear independence of the vectors (/Xi^o(O), . . . ,/ii/,r(0)) G W~^^, 
u = 1, . . . ,r, there exists A so that 

/loo(O) = . . . = /Io,m-l(0) = A, Jlom{0) = ■ ■ ■ = JioA^) = 0" 

Taking into account the condition mo(0) = 1 which is equivalent to 




(A(,o(0) + ... + //.,^-i(0)) = 1, 



we obtain A = 

10 



For the inductive step we assume that (10) is vahd for o; > and (a), (b) 
holds for all a' G Z'J^, a' < a. So, due to Proposition 6, there exist constants 

G C, 7 G Z'l, 7 < a such that (11) holds for all ^ < a. If 7 G Z'^, 7 < a, 
due to (5) and Leibniz formula, we have 



- y 



0</3<a-7 



a — 7 
/3 



m— 1 



^(27r^M-isfc)"-^-^DVfc(0) 



k=0 



x=0 



(15) 



It follows from (13) that 



y^Atofc/^t/fc = 0, = 1, . . . ,m - 1. 

fc=0 



Differentiating this equality a times gives 

r 

^D«-''Ai„t(0)DV*(0) = 0. 



E 

0</3<a 



fc=0 



Taking into account the inductive hypotheses, we have 



E 

0</3<a 



m— 1 



J] D"-/^/Iofc(0)DVfc(0) + ^>ofe(0)/i.fe(0) = 0. (16) 



fe=0 



k=m 



Multiply (15) by yjl^j and subtract from (16). After the same manipu- 
lation with each 7 G Z^, 7 < a, we obtain 

m— 1 

J]D-"/3^o/c(0)i^Vfc(0)- 



= 5^i^"/Iofe(0)/x,fc(0)+ J] 

k=m 

E („!J^ E 



0<7<a 



a — 7 



0</3<a-7 



0<^<a fc=0 

m — 1 



a — 7 

/5 



5](27rzM-is,)"-^-^L>Vik(0) 



fe=0 



J] L'"/Iofe(0)At.fe(0) + J] 

0</3<a 



k=m 



m— 1 



5^ (i?-/^m(o)- 



A;=0 



- E 



O<j<a-0 



a — 7 
/3 



q; — 7 / \/9 



q; 



A^(-27riM-isfe)°-/3-^ 



i^Vik(O). 
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Prom this, taking into account that 

a — 7\ / a \ f Oi\ ^ (a — py. f a — j3 



f3 J \a — ^ J \I3 J 7!(q; — /3 — 7)! \ 7 
and using the inductive hypotheses, the sum over (3 is deduced to 



(17) 



^ .. ^ / h—t\ 



0</3<a ^'^ ^ k=0 

-y—-^ ^ 

\ 7 / 

O<-f<a-0 / 



1 /« /5\ _^^_2^^^_i^^)„_^_, 

^ V 7 / 

0<7<a-/3 y 

fc=0 \ '^"^0<7<a^^^ / 



So, we have 



e(^>-«(o)-^e (;) 

k=G \ ^ 0<7<a ^ ' ^ 



A^(-27riM-isfc)"-^j /i^fc(O) + 

r 



k=m 



Similarly to the arguments for the initial step, it follows from (14) that there 
exists An, such that 



-7 — A; = 0, . . . , m — 1, 



m 

^"/iofe(O) = 0, k^m,...,r. 



Thus, (11) is valid ior P — a as was to be proved. 

Now we assume that (a), (b) are vahd. We will prove (10) by induction 
on a. If (11) is valid for a = 0, then /jtofe(O) = 1/y/m, k — 0, ... ,111 — 1. It 
follows from (13) and (b) that 

liuo{0) + ■ ■ ■ + llu,m-l{0) = 0, i/ = 1, . . . , r. 
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Hence, on the basis of (5), my(0) =0, u — 1, . . . , r, what proves the initial 
step. 

For the inductive step, we assume that (a), (6) is vahd for o; > and (10) 
holds for all a' G Zl, a' < a, i.e. 



D''-'^m^{M*~\) =0, 7eZ^, 7^0, 7 < a. 



x=0 



This yields (15) for 7 7^ 0. Multiply (15) by (^^"^ j and add to (5) 

differentiated a times. After the same manipulation with each 7 G Z^, 
7 < a, we obtain 



D'^m^{M*-^x) 

- y 



0<7<q; 



a — 7 



0</3<a-7 



0</3<a 



a — 7 



k=0 



P 



m—1 



J](27riM-is,)"-^-^L'Vfc(0) 



k=0 



- y 



0</3<a 



a 

a — 7 



a — 7 \ la 
/5 j V/J 



0<7<a-/3 

m—1 

J](27riM-^Sfe)"-^-^L'Vit(0). 

fe=0 



Due to (17), (11) and (6), this yields 



- E 



m—1 



\x=0 



0</3<a 



E E 

k=0 0<7<a-/3 



a — /j 
7 



E 

0</3<a 



-27r^M-isfc)"-/5-TDVfc(0) 

/ m—l 

^L'"-/^/Iofc(0)L'Vik(0) = i^" J]/Iofc(^)/^.ifc(i 



m—1 



fc=0 



.A;=0 



2:=0 



x] 



.fe=0 



a;=0 



It follows from (13) that D''m^{M*-'x) 



x=0 



as was to be proved. 
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It is not difficult to see that {V'jfe^} VM'^ property if and only if (10) 
is valid for all a E Z^, [a] < n. The following statement follows immediately 
fi:om Theorem 7 and Proposition 6. 

Theorem 8 Let n G Z_|_, r > m — 1, the functions i-ii,,k: l^vk £ 
v,k = 0,...,r, have derivatives up to order n at the origin, the matrixes 
A/" := {fJ'uk}l^k=o> ^ '■= {V-yk]l^k=Q satisfy (13); and let masks fho, mi, . . . ,m„_i 
be defined by (5). Then condition (10) is valid for all a e l/", [a] < n, if and 
only if 

(a) there exist A-y e C, 7 G Z'^, [7] < n, such that Aq = 1 and (11) holds for 
k — 0, . . . ,m — 1; 

(b) D''Jlok{0) =0, k^m,...,r for all^ e Z^, [7] < n. 

(n) 

Let n G Z+, we will denote by Loo the class of complex- valued functions 
which are in L^oiT'^) and have continuous derivatives up to order n at the 
origin. 

Lemma 9 Let jjiyQ, Ji^q G L^^ , i/ = 0, . . . , r, and 

r 

y^^i^uov^M = 1- (18) 

Then there exist functions fi^k, J^uk £ L'£\ ly = 0, . . . ,r, k = 1, . . . ,r, such 
that 



^(J^uifJ^uk = hi, k,l = Q,...,r. (19) 



i/=0 



Proof. Set 



r 



l^uo / \lT.\l^io\^, if WEI/^wPt^o, 
=0 V /=o 



p = 0, . . . ,r. 



1/^/m, if WEI/^/o|^ = 0, 



It is clear that the functions h'^^q are essentially bounded and 

EKor = i- (20) 



u=0 
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T 

It follows from (18) that E l/"/o(0)P ^ 0. So, /x'^q ^ = 0,...,r. 

Let us extend the unit vector /iQQ, . . . , ii'^q to a unitary matrix. Due to (20), 
there exist so that Ail,Qo(0) 7^ 1- We may assume that — (else we 
will interchange //'^^o //qq) extend this new vector to a unitary matrix 
and interchange its 0-th and VQ-ih. rows). Due to Householder transform, an 
extension to a unitary matrix may be realized by: 

/^ofc - /^fco . _— ' /^z/fc - % - _— , t^, fc - i,...,r. 

Because of (20), we have |/x^q| < \/l — l/^ooP, = l,...,r. This yields 
essential boundedness of the functions ji'^^^. Since 1 — /Xoo(O) 7^ 0, it follows 
that //'^^ e I/, /c = 0, . . . , r. Set 

J^uk ■■= l^'uk^ 1/ = 0, . . . , r, /c = 1, . . . , , r, 

Qfc (A*0A;) ■ ■ ■ ) A^rfc), /c = 0, . . . , r, 

Qo (a*oo, ■ ■ ■ , A*r-o), Qfe Qfe — QkQoQo, k — 1, . . . ,r. 

It is not difficult to see that the entries of Qk are in L^oo and QkQf = Ski, 
k,l = 0, . . . ,r. It remains to denote by fi^k the i^-th component of Qk-O 

Lemma 10 Let A be a class of complex-valued functions such that 

(i) if f, g & A, a, 6 e C then af -\-bg & A, 

(ii) if f,9 ^ A, then fg E A, 

and let A be a class of matrixes whose entries are in A. If any two n x 1 
matrixes Q, Q E A satisfying Q^Q = 1 can be extended to n x n matrixes 
J\f,J\f G A satisfying N"^N = In, then any two n x j matrixes M., M. E A, 
1 < j < n, satisfying M.^ M. = Ij. can be extended to n x n matrixes 
J\f,J\f e A satisfying M^N = In, 

Proof . We will prove by induction on j. The base for j — Ijs given. Let 
us check the inductive step j — 1 — > j. Let j xn matrixes M., M. E A satisfy 

Ai^A4 = Ij. Denote by Q^, the fc-th columns respectively of A^, Ai. Due 
to the statement of the base, the matrixes Qi, Qi E A can be extended to 



nxn matrixes N',N' E A satisfying W'^W ^ In, I'et Q'f., Q'f,, k ^ 2, ... , 
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denote the k-th columns respectively olN' ^N' G A. Fix a point x for which 
Q'i{x)^Q\{x) — 5ki. Since the vectors Q2{x), . . . , Q'ni^) form a basis for the 
orthogonal complement to Q'i{x) in R", we have 

n 

Qk(x) ^^aik{x)Q'i{x), k^2,...,j. 

1=2 

Similarly, 

n 

Qk(x) ^^aik{x)Q'i{x), k^2,...,j. 



It is clear that aik, aik G A and Yl ctikCtik' = ^kk', k,k' = 2, . . . ,j. Due to 

1=2 _ 

the inductive hypotheses, there exist functions C(ik,aik G A, I = 2,...,n, 
k = j + 1, . . . ,n, such that 

n 

'^aikaik' ^ 5kk', k,k' ^2,...,n. (21) 



1=2 



Set 



Qk ■=^oiikQ'i, Qk ■=^oiikQ'i, k = j + 1,. . . 



1=2 1=2 

Because of (21) and biorthogonality of the systems Qi, Qg, . . . , Q'^ and 
Qi, Q'2, Q'n, we obtain 

rp - — - 

Qi{x) Q,^{x) ^ Ski, k,l^l,...,n. 

To complete the proof it remains to introduce matrixes Af and J\f whose 
columns are respectively Qi, . . . ,Qn and Qi, . . . , Qn- 

Now we are ready to give a necessary condition for VM"' property. 

Theorem 11 Let dual wavelet systems {ipj'k}, {i^jk} generated by refin- 
able functions (p, (p whose Fourier transforms have derivatives up to order 

n at the origin^ and let the entries Uvk^J^vk of the corresponding polyphase 
matrixes be in L^oo . // {V^j^ } has VM"^ property, then there exist 

complex numbers 7 e Z'^, [7] < Aq = 1 such that (11) holds for 
k = 0, . . . ,m — 1. 
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Proof. Due to Lemmas 9, 10, the (r + 1) x m matrixes A1, A1 can be 
extended to (r + 1) x (r + 1) matrixes M,N such that their entries are in 

and M'^M = Ir+i- It remains to apply conditions (a) of Theorem 8.0 
Next let us discuss how to construct a dual wavelet system with VM" 
property generated by a pair of refinable functions ip, (p with masks mo-mo 
respectively. Due to Theorem 11, we know that the polyphase representatives 
of rho should satisfy (11). But this theorem said nothing about mo. Consider 
the following example. Let d = 1, M = m = 2, sq = 0, si = 1, /ioo = A^oo = 
Hqi = -ij, /xoi = (l — |sin27ra;). It is clear that (11) holds for n = 1 

(with Ao = 1, Ai = 0). Assume that there exists dual wavelet systems {i^^'^j}}, 



{^jk } with VM^ property (for {'ipji s) generated by these functions. This 
means that the polyphase matrixes whose first rows are respectively 

(/ioo, /^oi)) (fioo, Jj'Oi), satisfy (6) and such that (10) holds for the corresponding 
wavelet masks. Due to Lemmas 9, 10, the matrixes M.,M. can be extended 
to matrixes J\f,j\f such that their entries /i^k, V'vki i', k = 0, . . . ,r, are in L 
and N^^N" are mutually inverse. It follows from Theorem 8 that 



(2) 
oo 



But 



dx 




^/iofc(a;)/iofc( 



X) 



k=2 



0. 



x=0 



^l^ok(x)i^ok{x) 



x=0 



d f 11, 

— 1 1 sin 27ra; 

dx\ 2 2V 2 




1 -^l^ok{x)iJi,Qk{x 

k=0 

d f i 
= — - sm 2tix 
dx\A 



x=Q 




x=0 



(22) 



So, we see that a generating refinable function </? should be also chosen 
properly to provide VM"^ property for {ip^j'j}}- 

The following statement gives a sufficient condition for a pair of refinable 
functions to generate dual wavelet systems with VM^ property. 



Theorem 12 Let </?, (p he refinable functions, their Fourier transforms (p 
have derivatives up to order n at the origin, i^(0) = ^(0) = 1, and let 

jiQQ, . . . , jiQ^rn-i, A^OO) • • • ) Aio,m-i ^ -^c»'' he the polyphase representatives of 
their masks. If there exist complex numbers A^, 7 e Z^, [7] < n, Aq = 1, 
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such that (11) holds for k — 0, . . . ,m — 1 and there exist functions ^q^, 
Jiok G k — m, . . . ,r, such that 



fc=0 

D%k{0) = 0, k 



V/3 e Z^, [p] < n. 



then the functions (p,(p generate dual wavelet systems {ipj^}, {'ipjk'} with 
VM^ property for {ipj^}- 

Proof . Set Q = (yUoo, • • • , A^or), Q = (a*oo; ■ ■ ■ , J^or)- Due to Lemma 9, 
the 1 X (r + 1) matrixes Q, Q can be extended to (r + 1) x (r + 1) matrixes 
•A/" = {lJ'iyk}l,k=oi — {fii^k}l,k=o ^^'^^ ^^^^ their entries are in and 
AfAf^ — Ir+i- So, the matrixes 

M := 



( 



M 



satisfy (6). It follows from Theorem 8 that the corresponding wavelet masks 
mi, ... , rrim-i satisfy (10) for all a e if- ^ [a] < n, what was to be proved. 
Applied mathematicians and engineers are especially interested in con- 
struction of compactly supported wavelet systems. To provide this property 
generating refinable functions should be compactly supported and wavelet 
masks should be trigonometric polynomials. 

Theorem 13 Let (/?, Lp he compactly supported refinable functions with poly- 
nomial masks, (p{0) = ^p{0) = 1, and let yUoo, • • • , A^o,™-!; l^'OO, ■ ■ ■ , l^o,m-i 

be 

the polyphase representatives of their masks. If there exist complex numbers 
X^, 7 G Z'^, [7] < n, such that (11) holds for k = 0, ... ,m — 1 and there 
exist trigonometric polynomials /lok, J^ok, k — m, . . . ,r, such that 

r 
*;=0 



D^JiokiO) = 0, k 



V/? e zi, [p] < n. 



then the functions y?, ip generate dual compactly supported wavelet systems 
{^tk)' {'^jfc^} VW^ property for {ip^^^] . 



Proof. Set Q = (/ioo, ■ ■ ■ , A*or)) Q = (fioo, ■ ■ ■ ,J^or)- Due to Suslin's 
solution of a generalized Serre conjecture [18], the row Q can be extended to 
a unimodular matrix with polynomial entries. After this it is not difficult to 
find (r + 1) X (r + 1) matrices A/", A/" extending Q,Q. such that their entries 

are trigonometric polynomials and MM"^ = Ir+i (see [15], [8], [12, §2.6]). 
Next we repeat the arguments of the previous proof. 

Let us return to the example before Theorem 12. We could not succeed 



with VM^ property because the derivative of ^ lJ'Ok{x)lJ'Ok{x) did not vanish 

fc=0 

at the zero. Try to change /^oi = for /xqi = A*oi- Now we have 



d 

dx 



— - sm iTxx 
dx\S 



Kk=Q 



x=Q 



= 0. 

x=0 



So, providing condition (11) for both the masks mo, mo improved the situa- 
tion. Taking into account this observation, let us consider generating masks 
mo, mo whose polyphase representatives satisfy 



D^fiokiO) = ^ E ^7 (^] i-^niM-^kf-^ V/5 e Z^, [/5] < n, (23) 
D%k{0) = ^ E ^J^)i-'^^'^M-'skf-^ V/5 e [/3] < n. (24) 

Theorem 14 Let dual wavelet systems {'j/'j^''}, {i^jk} generated by refin- 
able functions (p, (p whose Fourier transforms have derivatives up to order 
n at the origin, the entries ^XukiV-vk of the corresponding polyphase matrixes 
be in and there exist complex numbers A«,, A^ , 7 G X'^, [7] < n, 
Aq = Aq = 1, such that (23), (24) are fulfilled for k = 0, . . . ,m — 1. If at 
least one of the systems {'ipj'^J}, {'^jk} VM"^ property, then 

E \] ^7^a-7 = Va e Z^, < [a] < n. (25) 



0<7<q; 



Proof . Let a e Z^^, < [a] < n. A; = 0, . . . , m - 1, p := 2'KM-^Sk. Due 
to (23), (24), we have 

mL*" (//ofe(x)^^) [^^ = E (")^''m(0):D^^^^^^ = 



0</3<a 
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E E E (;)(Q("/)(-irM-A,_,Ai;:;= 



E E E (;) (0 (" /) (-i^"")^"^^-^"-^- 

0<7<a: 0<(5<a-7 7</3<«-(5 / \ ' / \ ' 



^ E E (;)(Q(::^^)(-i)wv,w,.,= 



0<7<Q! 7<e<a 7</3<a— e+7 



EM-E E (;) (Q (::,^) 



0<e<a 0<7<e7</3<a-e+7 

a \/'/t + 7\/a — /« — 7 



EM'E E (.:jr; j .-7 )'-^'^'"^"-= 



0<e<a 0<7<e 0<ft<a-e 

a! 



0<e<a 0<K<a-e ^ ^ 0<7<e ^ ' ^ 

em^ e a.a:::z;( " e (^)^«^-(26) 

0<e<a 0<K<a-e ^ ^ Q<k<o. 



does not depend on k and 

1=0 



So, (^/xoik(a;)/xofe(x)) 

m— 1 

E (") A^Aa-K = E (/^oifc(2^)A«ofe(a:;)) = ( E l^ok{x)liok{x) 



0<K<a k=0 



' m— 1 



a;=0 



Due to Lemmas 9, 10, the (r + 1) x m matrixes Al, can be extended 
to (r + 1) X (r + 1) matrixes N^N such that their entries are in and 
H'^H — Ir+i- Due to conditions (b) of Theorem 8, 



= 0, k — m, . . . ,r. 

x=0 
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Prom this, taking into account that ^ l~iok{x)nok{x) = 1, we obtain 

A;=0 

m—l 

a;=0 



.A;=0 



Kk=0 



x=0 



Corollary 15 Let /xofc^, /xqa; G Dx! for some k — 0, . . . ,m — 1. If there exist 
complex numbers Xj,Xj, 7 e Z^, [7] < n, such that (23), (24), (25) are 
fulfilled, then 

-D" (m(x)ilok(x)) = 0, 

\ / x=0 

for all aeZ'l,0< [a] < n. 

The proof of Corollary 15 follows from (26). 

4. Methods for construction compactly supported wavelet 
frames with VM'^ property. 

Due to Corollary 15 and Theorem 13, refinable functions (p whose masks 
are trigonometric polynomials satisfying (23), (24), (25) generate dual com- 
pactly supported wavelet systems {i/^jk}, {V'jfe^} with VM'^ property for 

{V-'ifc } (we can set /xom = l,Aiom = Yl f^'oi f^oh this case). But this is 

_ 1=0 

a bad construction. The system {"ip^}} docs not have VM^ property. As was 
mentioned above, such a system can not be a frame. The following method 
allows to provide FM" property for each of the systems {'0^^'}, {i'j'k}- 

Step 1. Given n e Z*^ and given a set of parameters e C, /3 e Z^, 
[P] < Xq — 1, find a dual set of parameters A^ e C satisfying (25) by the 
following recursive formulas 



-^0 — 1; Xa — —Xa — {ft] ^P^"- 



0<f3<a 

Step 2. Chose functions //qq, . . . , //q^_^ and /loo, ■ ■ ■ , A*oto-i defined by 



[a]=n+l j=l 
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[a]=n+l j=l 

where T^^q, Tk^a ai'e arbitrary trigonometric polynomials, Qa are trigonometric 
polynomials such that £'"^a(0) = 1, £"^^a(0) = for all /3 G Z^, /5 ^ a, 
1(3] < n (recursive formulas for computing ga are given in [17]). It is clear 
that (24) are fulfilled 

m—l 

Step 3. Set (t := l^'oi J^oi, t^ok ■= (2 - (t)^^;,, k = 0, . . . ,m - 1. 

1=0 

Due to Corollary 15, we have D'^a{0) = for all /3 G Z^, < [p] < n. It 
follows that D/^nokiO) = £'^/Xofe(0) for all /? G Z^, < n. It is not difficult 
to see that (23) holds and 

m—l 



1 - E l^ok{x)lJ-ok = (1 - erf- 



k=0 



Set fj,om :— 1 — a, /lom :— 1 — cr. 
Step 4. Find matrixes 



M = 



\l^m,0 



A*0,m-1 A*0,m \ 
A*l,m-1 * 



M = 



fJ'0,m-l A*0,m\ 
A*l,m-1 * 



l^m,m—l 



J 



such that their entries are trigonometric polynomials and A^A^^ = /,„+!. 

Though the matrixes Ai.Ai can be constructed theoretically (see the 
proof of Theorem 13), it is very complicate to implement the algorithm in 
practice. Instead, we suggest the following explicit way (the payment of 
simplicity of this way is increasing of the redundancy) . 

Set //o,m+i = 0, Ato,m+i = 0. For each i/ = 1, . . . , m + 1, define 



A*i/fc '■— ^m+l-u,k "~ /^Ofc/^O.m+l-f) 



f^O,m+l-u, 

/J'uk '■= Sm+l~u,k ~ l^0kl^0,m+l-u, k = 0, ... ,771. 

"J+l .— [~ .\rn+l 



It is not difficult to see that the matrixes A4 := {A*i/A;}I7it^o' {/^i^fcllTfe^o 
satisfy MM'^ = Im+2- 
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Realization of suggested construction is simple, and there is a good chance 
to succeed in applications using the corresponding dual wavelet systems 
'ipj'ki'^jk- However, as it was mentioned above, we can be sure that these 
systems form dual frames only if additional conditions are fulfilled, in par- 
ticularly, if ifi, if decay fast enough. Checking of such conditions is very 
complicate in practice. Moreover, usually it is also difficiilt to check that 
ip,ip G L2(M'^). Situation is better for the case i/jj^,^ = ipjk- Namely, if 
(fi e L2(M'^) is a refinable function with a mask mo, lim^(^) = 1, periodic 

functions m,^, u — 1, . . . ,r, are so that the corresponding polyphase matrix 

^ A*oo ■ ■ ■ A*0,m-1 

y /^r,0 ■ ■ ■ /^r,m— 1 

satisfies 

M'^M^U (27) 

and = m^(M*~^a;)^(M*"^x), then {tjj'f,}} is a tight frame. This fact 

was presented in [16, Corollary 6.7] with an additional assumption which 
may be omitted if we improve the proof repeating one-dimensional arguments 
of [13]. So, for any appropriate polyphase matrix Ai, we may be sure that 
i'^jk} ^ tight frame whenever (p G L2(M'^) and ip is continuous at 0. If mo 
is a trigonometric polynomial, then (f is compactly supported, which yields 
continuity of (p whenever e L2(M.^). But this condition will be also fulfilled 
because, due to (27), we have 

m— 1 



E I'^o'^l' ^ 1' (28) 



fc=0 

and the following generalization of well known Mallat's theorem holds. 

Proposition 16 Let the polyphase representatives of a trigonometric poly- 
nomial mo satisfy (28), fiok{0) = l/y/m, = 0, . . . ,m — 1, and let f{x) : = 

oo 

n mo{M*-^x). Then f G L2{R'^), \\f\\ < 1. 

Proof. First of all we note (see, e.g., [12], §2.6) that 

m—l 

J2 \Mx + M*-'s)\'' = Yl \m{M*x)\'. (29) 

seD(M*) fe=o 
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It follows from (28) that 

\mo{x + M*-^s)\'^ <1. (30) 

seD{M*) 

k 

Set fk{x) = n ^o{M*~^x)xM*k[^i/2,i/2]d{^)- 

Taking into account that for any 1-periodic (on each variable) function g 

J g{x) dx — j g{x) dx 

[0,l]rf '■^■'^(^*)M*-i[0,l)d+M*-V 

(see, e.g., [12, § 2.2] ), we have 

Yl\mo(M*-^x)\'dx^m'' I Yl\mo{M*'-^x)\' dx ^ 

M*'=[-l/2,l/2]d [0,l]d 



m'' I \[\mo{M*^-'x)\^dx = 

„ k-1 

^ / \{\mQ{M*^~'x)\^\mo{x + M*-\)\^dx< 



-6-D(A^*)M*^i[0,l)d i=l 



^ fc-i „ fc-i 

k 



j Y[\mo{M*^-^x)\^dx^m^-^ j Y[\mQ{M*^-^-^ x)\^ dx 

-irn lid •?=^ rn.nd 



k-l 

\[\mo{M*-=x)\^dx=\\h. 

M*'=-l[0-l/2,l/2]<« 



l|2 
ill > 



Since, by (5), mo(0) = 1, the infinit product converges at each point, which 
yields fk{x) — > f{x) for all x e M*^. It follows from Fatou's lemma 

J \nx)fd^ = J fl\mo{M*-^0\'d^-supJ m)\'d^<^-0 

Now let us discuss how to choose a suitable refinable mask mo for con- 
structing compactly supported tight wavelet frames with VM"' property. 
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Note one more necessary condition for VM^ property. If {i/^j^} is a tight 
frame and Ai is its polyphase matrix with polynomial entries, then, due 
to Suslin's theorem[18] and Lemma 10, there exist (r + 1) x (r + 1) poly- 
nomial matrixes Af = {fi^kVu k=0' — {J^i^k}l,k=o extending M. such that 
NN'^ = Ir+i- If VM^ property holds for {ipj'k}^ W Theorems 8 and 14, 

D^^okiO) = ^ V A, (^] {-2mM-hkf-^ V/3 e Z^, [/?] < n, (31) 



0<7</3 

where a set of parameters A„ is so that 



Ao = 1, ^ [) A^A^-y = Va e Z^j., < [a] < n. 



(32) 



0<7<a 



Of course, providing (31) and (32) is not enough, in particular, these condi- 
tions do not guarantee (28). Now we will show that starting with arbitrary 
trigonometric polynomials satisfying (31), (32), it is possible to improve them 
such that the improved polynomials are polyphase representatives of a refin- 
able mask generating a tight frame with VM"' property. Our algorithm is 
based on the following result by M.A.Dritschel. 

Theorem 17 [4] IfT is a strictly positive trigonometric polynomial ofd vari- 
ables, d> 1, then there exist trigonometric polynomials ti, . . . ,tN so that 

T= + + (33) 



Note that another proof of this theorem was found by J. Jeronimo and 
M. J. Lai [5]. Their proof is constructive for the case d = 2. In the case 
d = 1, due to Riesz lemma, (33) holds with = 1 for any non- negative 
trigonometric polynomial T. 

Step 1. Given parameters A^, < [a] < n, which satisfy (32). choose 
functions //q^.. A; = 1, . . . , m — 1, defined by 

)"ofe(^) = ^ E 9a{x) f") (-27rzM-^Sfe))^A„_^ + 

[a]<2n 0<p<a 

J2 TUx)fl{l-e'-'^^r, ^ = l,...,m-l, (34) 

[a]=2n+l j=l 
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where T^^a ^-re arbitrary trigonometric polynomials, Qa are trigonometric 
polynomials such that D"g^{0) = 1, D^g^{Q) = for all /? e Z^, /? a, 
[p] < n (see [17]). 

Step 2. Choose large enough positive integers Li,Mi, i = 1, . . . ,d, such 
that the functions 



'M - n (1 - ^i"'^' ^^0"^' A*ofc(^)' k^O,...,m-l, 



i=l 



satisfy 



m— 1 



(35) 



k=0 



and L'^//(;fe(0) = D^fi',,{0), ^ = 1, . . . , d, for all /3 e Z^^, < n. 



ra—l 



Step 3. Set <j := E l/^o.^ 

fe=0 



Due to Corollary 15, L"^/xofc(0) = Df^^l^{Q) = D'^ix'^^^{0), k = 1, . . . ,d, ior all 
Step 4. Find trigonometric polynomials /lom, . . . , /iQr so that 



1 



fe=0 k=m 

D^fiokiO) = 0, A; = m, . . . , r, V/5 G Z^, < n. 
Such polynomials exist due to the identity 



(36) 
(37) 



1 



3 a 
2-2 ' ^ 



[1-a) 



£ > i 
4 — 2 

Step 5. Find a unitary matrix 



the fact that 1 — f > | (because of (35)), Theorem 17 and Corollary 15. 













M = 


A^io ■ 


■ A*l,m-1 


* 


* 




\ /^r,0 ■ 


fJ'r,m—l 


* 


■ * / 
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whose entries are trigonometric polynomials. 

It is not known if any appropriate row may be extended to a unitary ma- 
trix with polynomial entries. But if /ior = (which always can be realized), 
then M. can be defined by: for all = 1, . . . , r set 

fJ'vr '■— fJ'0,r-v, , fJ'uk '■— 5r-v,k — A*OfcA*0,r--i/, /c = 0, . . . , r. (38) 
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